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ABSTRACT 



Recent observations and theoretical work on gamma-ray bursts (GRBs) favor 
the central engine model of a Kerr black hole (BH) surrounded by a magnetized 
neutrino-dominated accretion flow (NDAF). The magnetic coupling between the 
BH and disk through a large-scale closed magnetic field exerts a torque on the 
disk, and transports the rotational energy from the BH to the disk. We investigate 
the properties of the NDAF with this magnetic torque. For a rapid spinning BH. 
the magnetic torque transfers enormous rotational energy from BH into the inner 
disk. There are two consequences: (i) the luminosity of neutrino annihilation is 
greatly augmented; (ii) the disk becomes thermally and viscously unstable in the 
inner region, and behaves S-Shape of the surface density versus accretion rate. 
It turns out that magnetically torqued NDAF can be invoked to interpret the 
variability of gamma-ray luminosity. In addition, we discuss the possibility of 
restarting the central engine to produce the X-ray flares with required energy. 

Subject headings: accretion, accretion disk- black hole physics - magnetic fields - 
gamma rays: bursts - neutrinos 
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1. INTRODUCTION 

The nature of the central engine of gamma-ray bursts (GRBs) remains unclear. 
Currently favored models invoke a binary merger or a collapse of compact objects. These 
models lead to the formation of a transient hot and dense accretion torus/disk around a 
black hole (BH) of a few solar masses. 

The typical mass accretion rates in GRB models are extremely high, of the order 
of a fraction of solar mass up to a few solar masses per second. Under such conditions, 
the disk becomes dense and hot enough in the inner regions to cool via neutrino losses. 
For this reason, Popham, Woosley & Fryer (1999, hereafter PWF99) named these disk 
neutrino-dominated accretion flows (NDAF). Energy extraction from the BH-accretion disk 
for powering GRB is also possible, such as by Blandford-Znajek (BZ) process (see Lee et al. 
2000 for a review of this model), and Blandford-Payne (BP) process and Parker instabilities 
in the disk (Narayan et al. 1992; Meszaros & Rees 1997). 

NDAF has been extensively discussed by many authors, e.g., PWF99, Narayan, Piran 
& Kumar (2001, hereafter NPK01), Kohri k Mineshige (2002), Di Matteo, Perna & Narayan 
(2002, hereafter DPN), Chen & Beloborodov (2006), Janiuk et al. (2004, 2007) and Gu, Liu 
& Lu (2006, hereafter GLL06). However, due to its low energy conversion efficiencies and 
the effects of neutrino opacity, the power produced by neutrino-anineutrino annihilation 
can hardly to match those of the energetic short-hard GRBs (e.g. GRB080913) and X-ray 
flares. Recently, some authors (e.g., DPN02, Fan et al. 2005, Shibata et al. 2006 and 2007 
and Perez-Ramirez et al. 2008) suggested that the MHD process should be considered in 
the disk model. Moreover, as shown in Shibata et al. (2006, 2007), the magnetic braking 
and magneto- rotational instability (MRI, Balbus & Hawley 1991) in the disk play a role in 
angular momentum transporting, which causes turbulent motion, resultant shock heating, 
and mass accretion onto the BH. On the other hand, researches showed that the magnetic 
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fields can be magnified up to 10 15 ~ 10 16 G by virtue of MRI or dynamo process (Pudritz 
& Fahlman 1982 and references therein) in hyperaccretion disk. These considerations 
stimulate us to discuss the magnetized NDAF. 

Based on the work of NPK01 and DPN02, Xie et al. (2007) discussed the BZ and 
BP processes in NDAF. They found that the jet of GRB may be magnetically-dominated, 
which is also obtained by MHD simulations of Mizuno et al.(2004). 

Recently, the magnetic coupling (MC) between the central spinning BH and their 
surrounding accretion disk has been paid much attention (e.g. Blandford 1999; van Putten 
1999; Li & Paczynski 2000; Li 2002; Wang et al. 2002). As a variant of the BZ process, the 
MC process exerts a torque on the disk, and transports the rotational energy from the BH 
to the disk. The effects of MC torque has been discussed in some disk models, for example, 
Lai (1998) and Lee (1999) in a neutron star with slim disk, Li (2002), Wang et al. (2002, 
2003), Kluzniak and Rappaport (2007) in a compact object with thin disk, Ye et al. (2007) 
and Ma et al. (2007) in a BH with advection-dominated accretion flow (ADAF). It is found, 
the disk properties are greatly changed and its luminosity is augmented significantly due to 
the rotational energy of BH extracted in the MC process. Therefore, it is attractive for us 
to investigate the effects of MC torque on NDAF. To highlight the effects of MC torque, we 
ignore other MHD process, such as BZ and BP mechanism, and we refer to this model as 
MCNDAF. 

This paper is organized as follows. In Sect. 2 we describe the MCNDAF model, 
which is a relativistic steady state thin disk. The effects of MHD stress are described by 
the dimensionless parameter a. The main equations are based on DPN02 and NPK01. 
Recently, GLL06, Chen & Beloborodov (2006) and Shibata et al. (2006, 2007) argued that 
the general relativistic (GR) effects are important for NDAF, so we introduce GR correction 
factors to the equations. The MC torque appears in the angular momentum equation. 
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We solve the set of equations for the solutions in MCNDAF in Sect. 3, and compute the 
neutrino and neutrino annihilation luminosities. Following GLL06, we include the neutrino 
radiation from the optically thick region in the computation for the neutrino luminosity. To 
show the effects of MC torque, we compare it with previous results. In Sect. 4 we discuss 
the stability of the accretion flow. We also discuss the physical origin of the instabilities in 
MCNDAF. Finally, we summarize our results and discuss some related issues in Sect. 5. 

2. NEUTRINO-DOMINATED ACCRETION FLOWS WITH MC EFFECTS 

Considering that central BHs are rapidly rotating in most candidate GRB engines, 
we discuss a model of a steady state disk around a Kerr BH, in which neutrino loss and 
transfer are taken into account. Our model is presented in the context given by DPN02, 
and the GR corrections are adopted from Riffert & Herold (1995). 

As mentioned by DPN02 and PWF99, although in GRB central engines the accretion 
rate may vary, it is expected to vary significantly only in the outer disk. Hence it seems 
reasonable to study the main properties in the inner neutrino-cooled disk by assuming a 
constant accretion rate. 

Because the gas cools efficiently, we are entitled to discuss the MCNDAF model in 
the context of a thin disk (Shakura & Sunyaev 1973). The accuracy of the thin-disk 
approximation is not perfect at large radii, where the disk is thick. Fortunately, the details 
of the outer region have little effect on the solution for the neutrino-cooled disk (Chen & 
Beloborodov 2006). 

The MCNDAF model is a relativistic steady thin disk, and the large-scale magnetic 
field contributing to the MC process and the small-scale tangled magnetic field related to 
the viscosity are included. We assume that these two kinds of fields work independently, 



- 6- 



and the large-scale magnetic field remains constant at the BH horizon. Following Blandford 
(1976) we assume that the magnetic field B z on the disk varies as B z ocC,~ n , where £ = r/r ms 
is the disk radius in terms of the marginally stable orbit r ms (Novikov & Thorne 1973), and 
n is the power law index indicating the degree of concentration of the magnetic field in the 
central region of the disk. 

Based on equipartition relation the magnetic field at the horizon is related to the mass 
density at the inner disk as follows (McKinney 2005), 



r>2 

-g^7 — P0,diskC , [i) 

where p 04isk = Mt g /r 3 , t g = GM/c 3 and r g = GM/c 2 . 

The MC torque is derived in Wang et al. (2002) based on an equivalent circuit given 
by Macdonald & Thorne (1982) as follows, 



™ in , / x f n/2 (I - P) sin 3 6d6 

T MC /T = 4a. 1 + q) / V [ ' . 2 (2) 

J 2-{l-q) sin 9 

where T w 3.26 x 10 45 ( T |f^) 2 (^) 3 # • cm 2 ■ s~ 2 , a, = Jc/(GM 2 ) is the dimensionless BH 



spin parameter defined by the BH mass M and angular momentum J, q = a/1 — a 2 and 
(3 = Vl D /VL H is the ratio of the angular velocity of the disk tt D = ((r 3 / GM) 1/2 + a^GM/c 3 )^ 1 
to that of the horizon tt H = a*c 3 /[2GM(l + q)]. 

The mapping relation between the angular coordinate 9 on the horizon and the radial 
coordinate £ on the disk is derived based on the conservation of magnetic flux as follows 
(Wang et al. 2003), 



cos#= / 6(a*;£,n)<i£ (3) 
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where f = (r/r^) 1 / 2 , 



8(a.;g,w) = - ^ 1 " ra ^V I+a *^" 2 + 2a *^" 3 ( 4 ) 

where Xms = ij- ms /r g ) 1 / 2 . 

A number of works (PWF99; Chen & Beloborodov, 2006) used accurate equations of 
relativistic hydrodynamics in Kerr spacetime to study NDAF. GLL06 found that the GR 
effect must be taken into account for the power of GRB. 

The relativistic correction factors for a thin accretion disk around a Kerr BH have been 
given by Riffert & Herold (1995), 



err c z r 



3GM ,GM ., ., 

5 = 1 -^- + 2 «*(^r) 3/2 ) (6) 



C=l-4^)3/ 2 + 3fe*, (7) 



c 2 r c 2 r 
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3xMc 2 
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4G 


s/rx 


(X 2 C A 


4 


V G 2 



alM 3 c 2 x ZalM 2 



D = I ^ ^ * — 8 dx. (8) 

The equation of the conservation of mass remains valid, while hydrostatic equilibrium 
in the vertical direction leads to a corrected expression for the half thickness of the disk 
(Riffert & Herold 1995; Reynoso, Romero & Sampayo 2006): 
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H^^I^M (9 ) 
w pGM\ C 1 ; 



The viscous shear t rtp is corrected by 



t rv = - a p^/A. do) 



The basic equations of MCNDAF are given as follows. 

1. The continuity equation: 

M = -2nrv r J:. (11) 

2. The total pressure consists of five terms, radiation pressure, gas pressure, degeneracy 
pressure, neutrino pressure and magnetic pressure: 

p = ii aT 4 + PkT { l+3Xnuc ) 27^^ /3^ /3 

12 m p V 4 ; 3 V 87rm/ >e 3 9 V ; 

where P ma5 = /3tP is the magnetic pressure contributed by the tangled magnetic field in the 
disk, and (3 t is the ratio of the magnetic pressure to the total pressure. u u is the neutrino 
energy density defined as (Popham & Narayan 1995) 

u v = (7/8)aT 4 Y W2 + 1A/3 

^t v J2 + 1/^3 + 1/(30 

In equation (fT3l) = T a , Vi + Ts^. is the sum of absorptive and scattering optical depths 
calculated for each neutrino flavor (z/ e , z/^, u T ). The absorptive optical depths for the three 
neutrino flavors are (Kohri et al. 2005) 
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r a , Ve ~ 2.5 x 10~ 7 T^H + 4.5 x 10- 7 T£X nW > 10 #, (14) 

T a ,u» = T a>UT ~ 2.5 X KT^if, (15) 

where X nuc is the mass fraction of free nucleons approximately given by (e.g., PWF99; Qian 
& Woosley 1996), 

X nuc ~ 34.8p- 3/4 T 1 9 / 8 exp(-0.61/T 11 ). (16) 
The total scattering optical depth is given by DPN02 as 

T S}Vi ~ 2.7 x lO-^pxoif. (17) 
3. Combining the conservation of the angular momentum with equation (ITT]) , we have 

-^-(MZ) + 47rr#MC = ^ = -^-(^r 2 t rv H), (18) 
cfr dr ar 

where I is the specific angular momentum of the accreting gas. The flux of angular 
momentum transferred magnetically from the BH to the disk, Hmc-, is related to the MC 
torque T MC by 



T MC = 4tt / H MC rdr. (19) 
Vanishing of t rtp (or g) at r ms leads to 
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GMD „ n „„ „ / A 



Mr 2 \l^- + T MC = g = -4irr 2 t rip H = 4nr 2 HaP\j (20) 
4. The equation for the energy balance is 



Q+ = Q~ (21) 

where Q + = Q V i S represents the viscous dissipation, and Q~ — Q u + Q p hoto + Qadv is the 
total cooling rate due to neutrino losses Q u , photodisintegration Q p hoto and advection Q a dv 
We employ a bridging formula for calculating Q u , which is valid in both the optically thin 
and thick cases. The expressions for Q u , Q p hoto and Q a dv are (DPN02; GLL06) 



o = y ( 7 /s^ 4 ) , 22) 

Wu ^(3/4)^/2 + 1/^+1/(3^))' 



Qphoto = 10 29 PioVr ^™ H erg ■ cm 2 s 1 , (23) 



<u = «AH a T 4 + |/*li±Z= + ^ )i (24) 

r 3 2 m p 4 3 

where 4u u /3 is the entropy density of neutrinos. Note that the cooling function given by 
the bridging formula reduces to the optically thin expression for small optical depths (as 
adopted in PWF99) but differs significantly from the latter at optical depths ~ I. 

By considering the MC effects, the heating rate Q vis is expressed as 



gQ! D _ 3GMM D T MC Q' D 
^ ms ~ Arrr Snr 3 B 4vrr ' { ' 



where the second term is the MC contribution. 
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As we can see from equation (1201) . the magnetic torque may deposit angular momentum 
in the inner disk, and this extra angular momentum must be transported outwards by the 
viscous torque in the disk, resulting in energy dissipation and increasing the disk luminosity 
based on equation (125]) . 

Defining Q G = 3GMMD/(8nr 3 B) and Q MC = -T MC Vt' D / {A-nr) as the contributions 
due to the gravitational release and the MC process, respectively, we have the ratio 
r] = Qmc/Qg versus the disk radius R = r/r g as shown in Figure 1. 

From Figure 1 we find that Qmc is much greater than Qg in the inner disk, where 
the neutrino cooling dominates. The ratio r] is very sensitive to the value of a* and n, it 
increases monotonically with the increasing a* and n. This implies that the MC effects are 
more important for the greater a* and n. For simplicity, we choose a* = 0.9 and n=3 in the 
calculations, and discuss the influence of their values in Sect. 5. 

We solve numerically equations (j!2p . (120]) and (T2TT) to find the disk temperature T and 
density p versus the disk radius with the given a*, n and m (where fn is the accretion rate 
in units of M s _1 ). We take X nuc = 1 for the fully photodisintegrated nuclear, which is 
appropriate in the inner disk. In the calculation, we do not include the cooling term arising 
from the photodisintegration, because it is much less than the neutrino cooling rate in the 
inner disk (Janiuk et al. 2004). Furthermore, a — 0.1, M — 7M and (3 t = 0.1 are adopted 
in calculations. 

3. EFFECTS OF THE MC TORQUE ON NEUTRINO ANNIHILATION 

LUMINOSITY 

As discussed in Sect. 2, the MC process applies a strong torque on the disk, resulting 
in huge viscous dissipation. This would lead to a more powerful neutrino radiation, and 
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neutrino annihilation luminosity. Here, we will show the effects of MC torque on the 
neutrino annihilation luminosity. 

To show this, we compare the results of MCNDAF with the NDAF model without MC 
(hereafter NDAF refers to the model without MC). GLL06 pointed out that the GR effects 
and the neutrino radiation from the optically thick region are important for the NDAF 
luminosity. Therefore we include these two effects in our calculations for both MCNDAF 
and NDAF. 

The neutrino luminosity from the accretion flow is 



We are interested primarily in the properties of the inner accretion flow, where neutrino 
processes are important. As argued in PWF99, NPK01 and DPN02, the flows are fully 
advection-dominated for r > 100r 9 , since neutrino cooling is not important and photons 
are completely trapped. Thus we concentrate the discussion in the region from r ms to 
?"max = 100r 9 . 

Our method for calculating neutrino annihilation is similar to PWF99 and Rosswog et 
al. (2003). The disk is modeled as a grid of cells in the equatorial plane. A cell k has its 
neutrino mean energy e k v% and luminosity I*., and the height above (or below) the disk is 
dk- The angle at which neutrinos from cell k encounter antineutrinos from another cell k' 
at that point is denoted as 6w Then the neutrino annihilation luminosity at that point is 
given by the summation over all pairs of cells, 




(26) 



k js-f k js- k j^i 




(l-cos6 kk ,) (27) 
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where A 1 ps 1.7 x 10 44 cm • ergs 2 • s 1 and A 2 ps 1.6 x 10 56 cm • erg's 2 s l . 

The total neutrino annihilation luminosity is obtained by integrating over the whole 
space outside the BH and the disk, 



As shown in Figure 2, the variations of L v and L vV versus fn for MCNDAF are indicated 
by the thin and thick solid lines, respectively, while those for NDAF are marked by the 
dotted and dashed lines, respectively. It is found that L u and L uD are greatly strengthened 
in MCNDAF. This means the spin energy acts as a powerful source for NDAF. 

According to our calculations for MCNDAF L vV varies from 3.7 x 10 49 ergs ■ s~ L to 
1.4 x 10 54 ergs ■ s -1 for 0.01 < fn < 10. We find that L uP nearly stays constant around 
~ 10 54 ergs ■ s" 1 for the accretion rate above fn ~ 0.5. This implies that the effect of 
neutrino optical depth becomes important. Our results for NDAF are in good agreement 
with those given by PWF99 and GLL06, but larger than those in DPN02. This is because 
the GR effects are taken into account in this paper as well as in PWF99 and GLL06. 

4. STABILITY ANALYSIS: THERMAL- VISCOUS INSTABILITY 

DPN02 discussed the thermal, viscous and gravitational stability properties of NDAF 
solutions. They found that NDAF is stable in most cases. But this result is not consistent 
with the variability in GRB lightcurve. To explain the X-ray flares, it is need that after 
the prompt gamma-ray emission has ceased, the central engine can be restarted (Fan & 
Wei 2005; Zhang et al. 2006). Based on this scenario, Perna et al. (2006) suggested that 
the X-ray flares could be produced by accretion of matter after the breaking of the disk 
due to the setting up of various instabilities either gravitational or viscous. Therefore, it is 




(28) 
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attractive for us to examine whether MCNDAF solution is stable. 
The condition for viscous stability is given by 

dM . . 

dE > < 29 > 

The stability curves for several radii in the disk are shown in Figure 3. 

From Figure 3, we find the rh — £ curves show S-Shape, in which the branch of 
solutions with negative slop is viscously unstable. It is shown that the viscous instability 
occurs at larger radius for larger accretion rate. However, we find this unstable can only 
occur at rh > 0.086, see Figure 5. 

It is clearly shown in Figure 3 that the disk is unstable at rh = 0.5 for R = 3, rh = 2 
for R = 10, and rh = 10 for R = 20. To understood this S-Shape we draw Figure 4. 

Inspecting Figures 3 and 4, we find the viscous instability occurs when the disk is 
neutrino cooing and radiation pressure dominated. For rh < 0.086, the MC torque will 
become very small, and the disk is optically thin to neutrinos. As discussed in NPK01, 
an optically thin NDAF is viscously stable for all pressure cases. Therefore, we will find 
no viscous instability when rh < 0.086. If the accretion rate beyond O.O86M s _1 , the MC 
torque transport enormous energy into the inner disk, and make the inner disk optically 
thick to neutrinos. In this time, we have rh oc X -1 for neutrino cooling and radiation 
pressure dominated case, and the disk will be viscously unstable. If gas pressure dominates, 
we have rh oc £ and mocS 3 for optically thin and thick NDAF, respectively. In the region 
where advection dominated, the disk is viscous stable, which is the well-known property of 
slim disk. 

From the above analysis, we conclude that the appearance of the viscous instability is 
due to the MC torque. First, the MC torque results in an optically thick neutrino-cooling 
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dominated flow in the inner disk part for high accretion rates. Second, the high opacity 
will drop the neutrino-cooling rate, and leave a hot and thick disk. In this region, the gas 
pressure drops and radiation pressure becomes important. We checked that the unstable 
solutions appear for black hole spin a* > 0.35 and magnetic pressure f3 t < 0.8. 

The disk is thermally unstable if (d\nQ + /dlnT) | s > (dlnQ~ /d\nT) | s . Then any 
small increase (decrease) in temperature leads to heating rate which is more (less) than the 
cooling rate, and as a consequence a further increase (decrease) of the temperate. For an 
optically thick NDAF in which neutrino cooling and radiation pressure dominates, we have 
Q + oc T 8 /S and Q~ oc T 4 . Therefore, the disk is also thermally unstable at the negative 
slope of the S-Shape m — E curves shown in Figure 3. 

It is noticed that Janiuk et al. (2007) obtain the viscous instability occurring at 
m > 10 without the S-Shape curves, which is caused by the behavior distribution for X nuc 
and photodisintegration term. However, in MCNDAF, the photodisintegration term is not 
included and X nuc = 1 is assumed for simplicity. Therefore, we infer that the S-Shape 
curves in MCNDAF arise from the MC effects. 

Finally, we check the gravitational stability condition, for which the Toomre 
parameter Q T should be larger than unity. For Keplerian disk, Q T is given by 
Q T = c s /t/(7rGS) = Q 2 D I '(nG p) . Q T decrease with increasing r so that the flow is most 
unstable on the outside. 

5. SUMMARY AND DISCUSSION 

In this paper we investigate some properties in MCNDAF. The angular momentum 
deposited in the disk by the magnetic torque exerted by the BH leads to a substantial 
additional dissipation of energy in the disk, which is greater than that expected from 
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gravitational release alone. Therefore, we obtain a series of MCNDAF solutions being 
different significantly from NDAF. 

The main results are summarized as follows. 

f . The neutrino annihilation luminosity in MCNDAF varies from 3.7 x 10 49 ergs ■ s -1 
to 1.4 x 10 54 ergs ■ s" 1 for 0.01 < rh < 10, while for NDAF the value range is from 
1.2 x 10 45 ergs • s" 1 to 2.6 x 10 53 ergs ■ s -1 , i.e., it is greatly improved by the MC torque. 

Recently, observations show that half of the Swift bursts exhibit X-ray flares. Fan et 
al. (2005) pointed out that the energy from NDAF cannot match the X-ray flares detected 
in GRB 050724 of ~100 s, which is also the time scale of the central engine. The time 
averaged isotropic luminosity of the X-ray flare component is Lx ~ 10 48 ergs ■ s" 1 . If 
we assume that the total mass available for accretion is ~ 1M (a typical value for the 
compact object merger scenarios and massive star collapse scenario), and that most of the 
mass is accreted during the X-ray flare phase, the time averaged accretion rate is about 
O.OlM s _1 . At this accretion rate, the jet luminosity powered by neutrino annihilation is 
L vD ~ 10 45 ergs ■ s -1 for NDAF without MC torque, and it is insufficient to power the X-ray 
flares. But the power produce by MCNDAF can satisfy this requirement. 

Very recently, the new observation of the highest redshift (z=6.7) swift source 
GRB080913 puts a very strong constraint on the central engine (Perez- Ramirez et al. 
2008). The duration of this short burst is T 90 = 8s, and the isotropic energy required is 
E iso 7 x 10 52 ergs. If the central engine is NDAF without MC torque, the jet collimation 
factor obeys fn < 7 x 10~ 4 , i.e., the jet should be strongly collimated. If invoke the MC 
torque in NDAF, the observed energy can be easily satisfied. 

2. The disk becomes thermally and viscously unstable in its inner region for m > 0.086. 
It is very interesting for us to obtain the m — £ curves behaving S-Shape, which may 
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produce a limit-cycle activity. The disk is rather thick in the inner region, and therefore the 
thermal and viscous timescales are close to each other. The timescale for the unstable can 
be estimated by the viscous timescale,^ = [1 / \aQ)](r / H) 2 , which is about 10ms at the 
inner disk. Following the discussions in Janiuk et al. (2007), these instabilities will lead to 
a variable energy output on millisecond timescales, which may correspond to the variability 
in the gamma-ray luminosity. The irregularity in the overall outflow can also help produce 
internal shocks. On the other hand, the thermal-viscous instability may be accompanied by 
the disk breaking, which can lead to the several episodic accretion events and explain the 
long-time activity accounting for the X-ray flares. 

Therefore, the MCNDAF can easily power both GRBs and its X-ray flares, naturally 
interpret the variability in the gamma-ray luminosity, and explain the production of X-ray 
flares. However, there are several issues should be addressed. 

First, in all of the MCNDAF solutions, we assume n=3 and a large BH spin a* = 0.9 
but does not give any reason. From Figure 1 we find the MC contribution is sensitive to 
a* and n. Thus for very small BH spin and value of n, the MC effects may be ignored, 
and the solutions return to NDAF solutions. Considering that the BH is spun-down in the 
MC process, while it is spun-up in the accretion process. The two processes with opposite 
effects result in a state with an equilibrium spin a% q . Calculation shows that a% q is greater 
than 0.85 for n > 3. This result implies that the MC effects are dominant in the whole 
duration of GRB, for which a fast-spinning BH is the central engine. 

Secondly, we made many simplifications in the MCNDAF model, such as we omit the 
photodisintegration term in the cooling rates, we assume the disk is thin, and a very simple 
magnetic configuration, and so on. Recently, Liu et al. (2007) took into account more 
realistic microphysics. Chen & Beloborodov (2006) worked out the NDAF solution under 
full Kerr metric. It is necessary to combine these effects with the MC process and work out 
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a more detailed model in the future. 

Finally, our MCNDAF is steady. Recently Janiuk et al. (2004, 2007) computed the 
time evolution of NDAF that proceeds during the burst. It is very interesting to investigate 
a time-dependent MCNDAF. 

Recently, Zhang & Dai (2007) proposed a hyperaccretion disk around a neutron star. 
They found, compared with a BH disk, the hyperaccretion disk around a neutron star can 
be cooled more efficiently and produce a much higher neutrino luminosity. As discussed by 
Kluzniak & Rappaport (2007), the magnetic dipole of the neutron star can also torque the 
disk. Therefore, it is also attractive to consider the effects of the magnetic torque in the 
context of hyperaccretion disk around a neutron star. 

We thank T. Liu for helpful discussions, and also thank the anonymous referee 
for his/her valuable comments and constructive suggestions". This work is supported 
by National Natural Science Foundation of China under Grants 10873005, 10847127 
and 10703002, the Research Fund for the Doctoral Program of Higher Education under 
grant 200804870050 and National Basic Research Program of China under Grant No. 
2009CB824800. 
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Fig. I. — The curves of the ratio of Que to Qg versus the disk radius R for a* = 0.9 with 
n = 3 (solid line), a* = 0.9 with n = 4 (dotted line) and a* = 0.8 with n = 3 (dashed line). 
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Fig. 2. — The total neutrino luminosity L v and L vD for a* = 0.9 and n — 3. The thick and 
thin solid lines represent L uD and L v of MCNDAF, respectively. The dashed line and dotted 
line represent L vV and L v of NDAF without MC, respectively. 
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Fig. 3. — The curves of m versus S for several given disk radii with R = 5,10 and 20 in 
solid, dashed and dotted lines, respectively. Other parameters are a* = 0.9 and n — 3. 
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Fig. 4. — Pressure (left) and cooling rate (right) components as a function of the disk 
radius, for three accretion rate values: rh = 0.5 (top), m = 2 (middle) and m — 10 (bottom). 
The pressure components are: gas pressure (solid line), radiation pressure (dashed line), 
degeneracy pressure (dotted line), and neutrino pressure (dot-dashed line). The cooling 
terms are: cooling rates due to neutrino emission (solid line) and advection (dashed line). 
Other parameters are a* = 0.9 and n — 3. 




Fig. 5. — The viscous instability is indicated by the shaded region in the parameter space. 
The parameters are a* = 0.9 and n — 3. 



